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Abstract. We determine the most general group of equivalence transforma- 
tions for a family of differential equations denned by an arbitrary vector field 
on a manifold. We also find all invariants and differential invariants for this 
group up to the second order. A result on the characterization of classes of 
these equations by the invariant functions is also given. 



In Lie theory, the invariance of functions and other objects under a transforma- 
tion group Q acting on an n-dimensional manifold V is usually characterized by 
the vanishing of functions under some vector fields generating the group action, 
and this vanishing is represented by a system of partial differential equations of the 
form 



where X — (x 1 , x n ) is a local coordinates system on V and where F is the 
unknown function. The importance of linear partial differential equations of the 
form (jl.ip usually referred to as determining equations for the invariant objects 
cannot be overstated. Indeed, they characterize invariant equations as well as their 
invariant solutions, and they have a similar importance in the study of Lie algebras 
and in representation theory. In physics, invariant operators of dynamical groups 
characterize specific properties of physical systems and provide mass formulaes and 
energy spectra [TJ [5] ■ Invariants of physical symmetry groups also provide quantum 
numbers useful in the classification of elementary particles [3! . It would therefore 
be desirable to consider the group of equivalence transformations of equations of 
the form and to determined all functions invariant under this group. Such 

functions are simply called invariants of the differential equation (jl.lj) . 

A method for the determination of invariants of linear and nonlinear equations 
build on an idea suggested by Lie himself [4] was developed in [5]- The method is 
based on the fact that the invariant functions for the infinite group of equivalence 
transformations of a given system of equations are precisely the invariants of the 
system of differential equations. That is, these functions are invariant under the 
group of transformations that confine the system of equations to a prescribed family 
of equations. The method also yields singular invariant equations, and has been 
used to complete the problem of determination of the Laplace invariants in [6] , and 
in [7] to characterize linearizable second order ode's. 



1. Introduction 




(1.1) 
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In the present paper we find the most general group G of equivalence transfor- 
mations leaving unchanged, except for its coefficients A%, an equation of the form 

n 

J2MX)d x ,U = (1.2) 

»=i 

where [x 1 , . . . , x n , U) S R" x M = M, and where we assume that none of the 
coefficients Ai for i = 1 , . . . , n vanishes identically. We then find the invariants and 
differential invariants up to the second order for this group and for an arbitrary 
number n of independent variables in the equation. We first treat with more details 
the cases n = 2, 3 before giving some generalizations of the results. Next, by 
investigating the regularity of the action of G on M, we show how the invariants 
found can be used to characterize families of equations of the form |L 

2. The group of equivalence transformations 



Owing to the linearity of equation (|1.2p . any invertible change of the dependent 
variable U and the independent variables (x , . . . ,x n ) = X that preserves the form 
of the equation should be of the form 

X = i/}(Y) (2.1a) 

U = H(Y)V(Y), H(Y) =£0 (2.1b) 

where Y = (y , . . . ,y n ) is the new set of independent variables, V is the new 
dependent variable and H is an arbitrary function. 

Theorem 1. The most general group G of equivalence transformations of equation 
(|1.2p consists of the set of all invertible changes of variables of the form 

x l = iP\Y) = ^(y% fori = l,...,n (2.2a) 
U = V. (2.2b) 

That is, each ip l {Y) is a function a the single variable y l , and G does not involve a 
change of the dependent variable. 

Proof. Under the general change of variables (|2. la|) . and by setting (j> = equa- 
tion (11.211 takes the form 



where 



Y,B J {Y)d y3 {U) = (2.3a) 

3 

Bj (Y) = X>(V00)§£wr)) = I>(*)§^P0 (2.3b) 

i i 

Equation (|2.3a.|) together with the expression of U given by (|2.1b[) shows that none of 
the coefficients Bj should vanish identically. However, the expression of Bj in (|2.3bj) 
shows that if ft (X) depends on more than one of the variables x 1 , for i = 1, . . . ,n 
it can be chosen as an invariant of an appropriate vector field, and so that Bj = 0. 
Hence (fp (X) = <fp(x q i), for some qj € {1, . . . ,n} and because of the invertibility 
of <p, <f> 3 must be a nonconstant map and all the variables x q3 must be distinct for 
j = 1, . . . , n. This implies in particular that x l = ^(y^) must also be a nonconstant 
function of a single variable. If we let a be the permutation that maps the ordered 
set {y 1 , . . . ,y n } onto the ordered set {y kl , ■ ■ ■ ,y kn } , then the «th component of 
ip o a~ l depends exactly on y 1 alone. On account of the arbitrariness of ip, we may 
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replace ip, by ip o a , and thus that we may always assume that x l = %p l (y l ), and 
equivalently y % = <p l (x l ). This reduces the expression for Bj(Y) in ([2~3b1 to the 
form 



where ip jl = dip j /dy j . Substituting (|2~Tbf into (|2.3aP and expanding, equation 
(11.21) takes the form 



J2H B S d yi V + vi^Bj d yJ h\=0. (2.5) 



The fact that the coefficient of V appearing in (12. 5| must identically vanish and 
the arbitrariness of the n coefficients Aj in the expression of Bj in (|2.4j) show that 
dyi H(Y) = 0, for all j = 1, . . . ,n. Thus H{Y) ^ is a constant function and 
without loss of generality we may assume that H = 1 . This last equality transforms 
equation (|2.5[) to the form 

J2B j (Y)d yi V(Y)=Q, (2.6) 

3 

which is of the prescribed form. This completes the proof of the theorem. □ 

Remark . It should also be noted that under the general change of variables (|2.1ll , 
it is always possible, by the well-known result on the rectification of vector fields, 
to put (|1.2[) in the form 

d y i(HV)=0, that is, (d y iH)V + H{d y iV) = 

Thus if we allow some of the coefficients Ai to vanish, then the only additional 
condition to be imposed on the change of variables (|2.1|) would be d y iH = 0, and 
all equations of the form (II. 2|) would be equivalent. There are clearly no invariant 
functions or invariant equations of any order in such case. 

We now move on to determine the infinitesimal generators of the group G. As 
already noted, equation (|2.2a[) implies that y l = </> l (x 4 ), for i — 1, . . . ,n, and this 
shows that the infinitesimal transformation of (12.211 has the form 



y^x' + eCix'), VkU, (2.7) 

where the functions £ l are also arbitrary, due to the arbitrariness of the functions 
ip 1 . The first prolongation of this transformation has the form 

d y iVnd x iU + e(-e'd x iU), (2.8) 

which implies that 

d xi U*d y iV + e(e'd yi V,) (2.9) 

where d x is the differential operator d/dx, for any variable x. A substitution of 
equation (|2.9[) into the original equation (|I.2p yields the infinitesimal transformation 
of that equation in the form 



E 



(Ai + eAiC ')d y iV = 0. (2.10) 
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This shows that the infinitesimal transformation A of the coefficient Ai is given 

Ai ~ Ai + tAi^ 1 . 

The infinitesimal generators of the equivalence transformation G therefore has the 
form 



3. Zeroth-order invariants 

We would like to first recall very briefly certain elementary facts about the invari- 
ant functions of a given transformation group. Suppose that the infinitesimal gen- 
erators of an r-parameters group of transformations G acting on the Q-dimcnsional 
manifold M are of the form 

Vl y,C J <>, ■ farfc=l,...,r. (3.1) 

3 

The invariant functions and invariant equations of G are determined by 

V k (F) = (3.2a) 

V fe (F) = (3.2b) 

F =0 

respectively, for k = 1, . . . , r. The number of fundamental invariants of G does not 
exceed Q — r, where r is the rank of the matrix J ') fc . of coefficients of the r 
operators Vk- Each of these functions naturally gives rise to an invariant equation. 
Invariant equations F = 0, where F is not an invariant function, and obtained 
by imposing the additional condition r < Q to the second equation of (|3.2b[) are 
often referred to as singular invariant equations. Using a Lie linearization test, 
such equations were recently shown [7] to characterize all linearizable second order 
ordinary differential equations. 

When some of the independent variables x J in the expression of the Vk can 
be taken as depend variables for other objects such as a differential equation, the 
generators Vk can be extended to involve higher order derivatives of the dependent 
variables. If V is a given infinitesimal generator of G, then we shall often use the 
same symbol V to represent both V and its m-th prolongation y( m ). Similarly, the 
mth jet space of M will often be denoted simply by M. 

Since the general change of variables (|2.2j) is merely a change of the independent 
variables and does not involve the dependent variable U, this variable is trivially 
an invariant for G. We shall therefore ignore this variable in our search for the 
invariant functions of G whose general form for the zeroth-order operator (|2.1ip is 
F(x\...,x n ,A 1 ,...,A n ). 

Theorem 2. The group of equivalence transformations G of (|1.2jl has neither 
invariant functions nor invariant equations. 

Proof. Rewriting the generic generator V in (|2.1ip as a linear combination of the 
arbitrary functions £ l and their derivatives gives 
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and this proves the first part of the theorem at once, on account of the arbitrariness 
of the functions £ l . To show that G has no invariant equation, we use an elementary 
technique similar to that used in [7j. Suppose that F(x , . . . , x n , A±, . . . , A n ) = 
is a nontrivial invariant equation for G, and so it explicitly involves at least 
one of the variables, say x 1 , in the set {a; 1 , . . . , x n , Ax, . . . , A n } . Then solving the 
equation for x 1 reduces it to the equivalent form x l — K(x 2 , . . . ,x n ,Ai, . . . , A n ). 
The arbitrariness of the functions and their derivatives implies again that we 
must have in particular 



cV (x 1 - K) 



= 0. 

c 1 =K 



But this last condition cannot hold because d x i (x 1 — K^j = 1, and this completes 
the proof of the theorem. □ 

4. First-order differential invariants 
The first prolongation V^ 1 ^ of the infinitesimal generator (|2.11[) of G has the form 

V" = V + E + EE- 4 >< <« J ' " ?') ^" (4-la) 

i i j^i J 

where we have used the notation 

OA. 

In terms of the linear combination of the arbitrary functions and their derivatives 
this expression takes the form 



A v = UT7> for i,j G {!,..., n}. (4.1b) 



n a n o 

^ dx l ^ OA 

I 

n 



(4.2) 



Equation (|4.2p clearly shows that any first-order differential invariant of G should 
be independent of all the independent variables x l , as well as the variables An . This 
last condition reduces V*- 1 -* to the form 

n 

V^=J2e'V^ (4.3a) 

i 

where 

V^=A^ + E(^-^)- (^b) 

As we are considering the arbitrary functions Ai in (|1.2p in their most general 
form, we may assume that the functions Aij — d x jAi do not vanish identically. It 
then follows from f|4. 3|) that the generators of the first prolongation of G depend on 
Q = n 2 independent variables. If however we assume that exactly p of the functions 
A^ vanish identically, then this number p is invariant under the action of G and 
Q = n 2 — p. 



Theorem 3. Consider the n operators V^n of (|4.3 
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(a) The rank of the coefficients matrix Ai of the operator V^n is n, which is 
maximal. 

(b) The V^n form an n- dimensional commutative Lie algebra. 

(c) The number fundamental first-order differential invariants of the group G 
of equivalence transformations of equation (|1.2p is n(n — 1). 

Proof. In any coordinate system of the form {Ai, . . . , A n , . . . } on the extended jet 
space on which the first prolongation of G operates, equation 1)4. 3ft shows that the 
first n columns of Ai is represented by the matrix diag {A\ , . . . , A n } which has rank 
?i, owing to the fact that none of the coefficients Ai is zero, and this proves part 
(a) and shows that the n vectors Ve*> are linearly independent. For part (b), if for 
any k G {1, . . . , n} we write 

V^ = A k d Ak+ ^(A kq ^-A qk ^-) 

then we readily see that the commutator [Vj>/, is a linear combination of the 
identically vanishing commutators 

[AijdAi, - ■[,.<> i i A k q dA kq - A qk d Aqk ] , [Aid Ai , A k d Ak ] 

[Aid Ai , A kq d Akq - A qk d Aqk ] , [Aijd Aij - .\,,t) \ , A k d Ak ] , 

This fact together with part (a) proves (b). Since the number of independent 
variables involved in the complete system of n operators Ven is n 2 , the number 
of their functionally independent invariants is precisely n 2 — rank(A^), which is 
n(n — 1). □ 

The most practical way to find the n 2 — n first order differential invariants of G 
would be to compute these invariants for low dimensions of M, i.e. for n = 2, 3 and 
then make use of the symmetry inherent in equation (|1.2p to find the invariants in 
the general case. 

For n — 2 and n — 3 we write equation p.2[) in the form 

aU x + bU y = 0, and aU x + bU y + cU z = 0, (4.4) 

respectively. In case n = 2, the operators Vga are given by 

V ? i / = ad a + a y da v - b x db x , , = bd b - a y d ay + b x db x ■ 

Solving the system of equations V^n(F) = 0, for i = 1, 2 by the method of charac- 
teristics shows that G has a fundamental system of invariants consisting of the two 
functions 

m Q*v b „ b T a 

T 12 - and Tai = 

a b 

In case n = 3, the three operators V^n are given by 

V ? i/ = ad a + (a y d ay - b x d bx ) + (a z d a ^ - c x d c J 
V ( 2, =bd b + (b x d bx - a y d ay ) + (b z d hz - c y d Cy ) 
= cd c + {c x d Cv - a z d a J + (c y d Cy - b z d bz ) 
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and the corresponding set of six invariants is found to be 



Tl2 = 


Cly b 

a 


T21 = 


b x a 


T13 = 


a z c 
a 


T31 = 


c x a 
c 


T23 = 


b z c 
b ' 


T32 = 


Cyb 
c 



The form of the invariants found for n — 2,3 together with part (c) of Theorem [3] 
asserting that the number of invariants in the general case is n(n — 1), which is 
2 (2) , suggest that all invariants can be found by associating with each subset of two 
elements of the set of n coefficients of the differential equation a pair of invariants 
according to a very simple rule. 

Theorem 4. The n(n— 1) fundamental invariants of the group G of equivalence 
transformations of equation (|1.2ll are given by 

Tii = ^, fori ™here .1,, (4.5) 

and where Ai and Aj run over the set coefficients of the equation. 

Proof. It is easily verified that the identity V^n(Tk g ) — holds for all i = 1, . ..,n 
and for all k 7^ q. Next, the functions Aij for i,j S {l,...,n} are functionally 
independent by assumption, and each depends on exactly one of them. □ 

Note that if we restrict the action of G to a sub-family of equations of the form 
(II. 2p for which exactly p of the function A^ vanish identically, then the maximal 
number of functionally independent first order differential invariants is n(n — 1) —p. 

5. Second-order differential invariants 
The second prolongation of the generator (|2.1ip of G has the form 

3 

+ [(?'-2?')A jii -A ji ?"] 8 Ajii 

+ 2 £ (?'-?' -^') A jik9 Ajik , (5.1a) 



i<k 



where as usual 



_ d Aj d Aj 

<' 77^' Ajlk = d^d^' etc - (5 - lbj 

Rewriting this expression as a linear combination of the arbitrary functions and 
their derivatives shows that any invariant function should be independent from the 
independent variables and from variables of the form Am for i = 1, . . . ,n. This 
reduces the expression of V*- 2 - 1 to the form 

n 

V (2) = j^e %i> +e (5.2a) 
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where 

V e , — Aid Ai + X An - A H d A + a vo9a w ^A^Ma J 

V k / 

+ 2 ^ A yfe ^ ijfc (5.2b) 
J<fe 

V fi „ = Ai 8 Au + ( 2A ij 9 Auj - Aji d AjH ) . (5.2c) 



It readily follows from equations (|5.2| that the second order differential invariants 
of G depend in general on n+n( n ^ 2 ) —2n variables, that is on n 2 (3 + n)/2 variables. 
This is the dimension of the subspace of the extended jet space of M on which 
the second prologation of G acts. 

Theorem 5. The set of operators {V^*/}™^ and {Vfi"}._j given in (|5.2p each 
generate an n dimensional commutative Lie algebra. 

Proof. Thanks to the term Ai&Ai appearing in the expression of each generator Vj; i 
as the only term involving d At , the coefficients matrix of these operators admits a 
submatrix of the form diag {Ai, . . . , A n } , which is clearly of rank n, showing that 
the V^i' generate an n-dimensional space. Similarly, as the term Aid An appears in 
the same manner in the expression of each generator V^m, the set {Vf*"}'^, also 
generate an n dimensional space. For each pair {i, k} , it is easy to see as in the 
proof of Theorem [3] that each of the commutators [V^>, V^ki] and [V^i»,Vtk«] is a 
linear combination of identically vanishing commutators. This completes the proof 
of the theorem. □ 

If we denote by £*w the jth derivative of £*, then Theorem [5] asserts that for 
j fixed, the Vj;(j)'s form a commutative Lie algebra for j = 1,2. However, the set 
of all operator V^u) for i = 1, . . . , n and j — 1,2 that determine the second order 
differential invariants of G do not form a Lie algebra in general when they are 
considered together, as this easily appears from the low dimensional cases. 

Indeed, if for n = 2, 3 we rewrite equation (|1.2|l as in (|4.4[) . then for n = 2, we 
have 

d d d d d d 

= a ?n + a yjrr + a yy~xZ, bx 7rr ~ 2bxx ~ai — 2bx y 



da da y da yy db x db xx db xy 

d d „ d d , d 



V « 2 ' = b db a *da- y ~ 2axy daV y ~ 2a ***°»» + b *db~ x + bxx db. 
d d d 

= a a 1" 2a V a oT — 

da x da xy db xx 
V ?2 // = -a y - h b— h Zb x 



d a yy db y d b xy 

In this case we have 

[V 5 i V e »] = V e „ , and [V f 2 , , V 5 2 „] = V e ,, 

However, the span of {Vgi/, V^>i V^i», Vja//} does not contain the commutator 
[V^u), V^k(p)} for any sets {i,k} and {j,p} of distinct elements. For instance, we 
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have 

9 

[Vji (2) , V^2 (i)] = — 2a y — . 

O (X x y 

We have a similar situation in the case of three independent variables. The opera- 
tors V^i (3) are given in this case by 

9 9 9 9 „ 9 9 



Vfi' = a- h a y - h a z - h a yy - h 2a yz - h a zz - — 

da da y da z da yy da yz da z 

9 9 9 9 9 

VxTTj £u xx — ^VxyTTi ^b xz — Cx ~ 

db x db xx db xy db xz dc x 



9 9 9 

— 2c xy — 2c ^qT- 

u i^xx u L -xy u ^xz 



,0 ^r,O d d d 

^? 2 ' = "TT - a y^ 2a xy - 2a aa - 2a yz - h b x —— 

db da y da xy da yy da yz db x 

9 9 9 9 9 

"T b z — ~t- b xx — -\- 2b xz — -t- b zz — Cj/n 
db z db xx db xz db zz dc y 

_9 _9 9 



xy q yy fi y z c\ 

o c xy o c yy o c yz 

9 d d d d.d 

dc da z da xz da yz da zz db z 

9 9 9 9 9 

~ db yz db db~ + Cx fhT + Cy d~^ 

9 _9_ 9 



xx a ^ x v pi ' yy fx 
o c X x Q Cxy (J Cyy 



d d d d d 
Vfi" = a— 1- 2a v — h 2a 2 — b x — c ; 



i~i 1 ^ '~* j y <~\ i*- ili 'Zi~i '-'x i~i 

^a;a; C C xx 

V ^+b^+2b^+2b^ ^ 

^ v da yy db y x db xy 2 db yz V dc yy 

9,9 9 9 9 

V^3 // = — a 2 ~ b z — h C- h 2Ca: ~ 1- 2c v 



oa zz db zz 0c z dc xz dc yz 

As in case n = 2, we have [V^i > , V ? i »] = " . That is, { V^i < , // } spans a solvable 
Lie algebra with nilradical {V^m } , for i = 1, 2, 3. However, here again the span of 
{V^i(j)} .^. does not contain the commutator [Vei y) , V^k ( P )] for any sets {i,fc} and 
of distinct elements. Indeed, we have for instance 

9 

[V^i (2) , Vj3 (i) ] = 2c x — . 

O Cxz 

There is no guarantee in this case that the number of invariant attains its maximum 
which is Q — t, with the usual notation. More over, they are much more difficult 
to find using the method of characteristic, we shall therefore attempt to determine 
the invariants of the second prolongation of G using the so-called method of total 
derivatives [HUH]- 

Suppose that we are given a system of equations of the form (|3.2ap where the 
Vfc's are arbitrary linear differential operators given as in (|3 . 1 [) and depend on a 
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total of Q variables. Denote again by r the rank of the coefficients matrix (£ ■ , 
and set p = Q — r. Thus we can solve (I3.2al) for r of the variables d Xt F in terms of 
the remaining p others, and this gives rise to the Jacobian system 

A t F = ^ + J2Us,t^=0, fort = l,...,T, (5.3) 

s—1 

where we have renamed the remaining p variables x T+ j as Uj, for j = 1, . . . ,p and 
where the /7 Sj t's are functions depending in general on the Q vriables X\, . . . , x r and 
Mi, . . . , Up. In this case the equivalent adjoint system of total differential equations 
takes the form 

r 

du s = *^2u s , t dx t , for s = 1, . . . ,p. (5.4) 

t=i 

The equations (13.2a|) and (|5.4j) are equivalent in the sense that they have the same 
integrals [8]. We denote by M. u the coefficients matrix {£/,,*} that completely 
determines the adjoint system (|5.4[) . 

For n = 2, by permuting the coordinate system on M so as to have in diag {a, 6, a, 6} 
as the submatrix corresponding to the first 4 columns of the coefficients matrix for 
the system S22 = { Vjri / , Vgi 1 , V^i » , } , we obtain the transposed matrix M. ^ of 
M. u in the form 



Ml = 



The corresponding system (15. 4p of total differential equations can be solved using 
methods described in [8] . We get stuck with a problem of finding some integrating 
factors while trying to solve by the method of characteristics the equivalent system 
(|3.2ap of linear partial differential equations for the system of operators S22 ■ How- 
ever, we readily get the following set of six functions by solving the corresponding 
adjoint system (|5.4[) . 



/ 





a yy 
a 


a 


a y 
a 


_2b^ 
a 


2^y\ 

a 




2a xy 


2a vy 


b., 


a y 









b 


b 


b 


b 


b 




2a y 


















a 








a 




V 





a v 
b 











2b x , 
b ' 



J 12 — , ^12 — \-o yi J\2 = £a x 

a a v a 

b x a o xx a o xy ab 

J21 = —r- j K 2i = — 1 1- ii, -J21 — — ; ^o y . 

b b x b x 

Although their number corresponds to the maximal number of functionally inde- 
pendent invariants in this case, not all of them are actually invariants because the 
system S22 is not complete. More precisely, only TJy and Kij, for i,j — 1,2 are 
invariants, and not only the J^'s are not invariants, but also the equations Jij = 
are not invariant equations. 

Similarly for n — 3, the total number Q of variables defining the invariants is 27, 
and we have r = 6. By permuting again the coordinate system on M, so as to have 
diag {a, 6, c, a, 6, c} as the first six columns of the coefficients matrix for the system 
of perators £32 = { V^i 1 , V^2 , , / , V^i » , V^2 » , V^3 » } , we obtain a more convenient 
representation of the 21 x 6 matrix M. u . The transpose M. 1 ^ of M. u in which only 
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its first six columns are represented has the form 
/ 



mT = 



a yy 2a yz a Z z _ ^ \ 

n n n n . ' 



\ 









a yy 


2a yz 




_b 3L 


a 


a 


a 


a 


2a xy 





2a yy 
b 


2a yz 
b 





l>.r 
b 





2a xz 
c 





2a yz 
c 


2a zz 

c 





2a y 
a 


2a z 
a 




















a y 
b 





























where the dots represent the remaining 15 matrix columns. Solving the correspond- 
ing system (|5.4|) yields the expected maximal number of 21 functionally independent 
functions. But since here again the corresponding system of operators 1S32 is not 
complete, only 15 of them are actually invariants of G. Reverting back to the orig- 
inal notation Ai = a, A2 = b and A3 = c, and using (|5.1b[) . these 15 invariants can 
be written in the form 



T ■ — 



A,, 



ij ~ A - 
-fin 



A 



j j ■ 



Lijk — A 



'<]!< 



AjAk 
~A~ 



(5.6) 



where i,j 6 {1,2,3}, with i 7^ j, and where {j, k} = {1,2,3} \ {i} for i = 1,2,3. 
We have thus obtained the following result. 

Theorem 6. Let TV be the maximal number of functionally independent invariants 
of the second prolongation of the group of equivalence transformations of (|1.2|) in 
n independent variables. 

(a) For n = 2, J\f = 4, and the invariants are the function Tij and Kij of (|5 . 5[> . 

(b) For n = 3, Af = 15, and the invariants are the functions Ty, Kij and L^k 
given by (|5.6|) . 

Contrary to the case of the first prolongation of G, a determination of all in- 
variants of the second prolongation for larger values of n using only invariants of a 
lower order of n and a symmetry argument does not seem to be obvious. Indeed, 



the equations (|5.5p and (|5.6p show that for n = 3, the invariants of type and 
can be simply derived by symmetry from those for n = 2 without any further 



calculations. However, the invariants of type L^k in (|5.6|) cannot be obtained from 
(|5.5p using only a symmetry argument. This makes it more difficult to find all the 
invariants for the second prolongation of G when n > 4. Nevertheless, we do have 
the following result which is solely based on a symmetry argument. 



Theorem 7. For n > 3, a fundamental set of invariants of the second prolongation 

Jijk of (|5.6p . whose total 



of G includes all the invariants of type Tij, Kij and Li 



number is n(n + n — 2)/2. 

Indeed, this result clearly follows from (15. 6p and the symmetry inherent in (II. 2p . 
by noting that the total number of the Tij, Kij and Lijk f° r ^ > 3 is 



+ 3 



i n(n 2 + n-2). 



(5.7) 



Although we may not find all the invariants of the second prolongation of G for 
larger values of n > 3 using only symmetry arguments, it should be possible to 
predict their number. Denote by M^' J the number of fundamental invariants of the 
fcth prolongation of equation Q1.2p (with n independent variables) involving terms of 
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the form Ai, where I is an index of the form i\t2 ■ ■ - ij with distinct i& £ {1, . . . , n} 
for k = 1, . . . , j. Note that the corresponding type of functions appears for the first 
time as invariants of (|1.2p when the number of independent variables is j, where 
2 < j < n. If we also denote by MJ} the number of invariants of the fcth prolongation 
for n variables, then a closer look at equation (|5.6[) suggests that A\ J = j(") and 
= M™ + W n , where 

w n = + + . . . + = ( n ) ■ 

j=2 ^ J ' 

Using the properties of binomial coefficients, it can be shown that ^2 j ec l ua l s 
n(2™ -1 — 1). Since by Theorem[3]we have M™ = n(n — 1), our conjecture follows. 

Conjecture. For any value n of independent variables in equation (|1.2p . the num- 
ber M£ of functionally independent invariants of the second prolongation of G is 
n(2"- 1 +n-2). 

This conjecture says that M| = 40, and Mf = 95. By a result of Lie (see [5]), it 
is possible to find differential invariants of G of higher order than 2 using invariant 
differentiation, but we will not discuss that here. 

6. Properties of the invariants 

It follows from Theorem[T]that every element of the equivalence transformations 
group G of equation (|I.2p can be represented by an invertible map <p of the form 

0: R n -> R n : X = (x 1 , . . . , x n ) ^ Y = <j>{X) = {^{x 1 ), (j) n {x n )). 

That is, the ith component of 4> depends only on the single variable a; 1 . In a given 
coordinate system X — (x 1 , . . . , x n ), each equation of the form (|1.2|) can be rep- 
resented by the n-tuple (Ai(X))™ =1 or just (Ai(X)). It follows from equation (|2.4p 
that under the action of 4> S G, the differential equation (Ai(X)) is mapped to the 
differential equation (Bi(Y)) = <j> ■ (Aj(X)), where 

Bl (Y) = MrHYWmy)) = 

and where ip = <f>~ x ■ If 9 is any other element of G, and Id is the identity transfor- 
mation, it is easy to see that 

Id-(^pO) = (MX)), and 6 ■ (<j>-(A i {X)))=6o<j>.(A i (X)). (6.1) 

Thus if we denote by E n the variety of all differential equations of the form (| 1 .2(1 , 
then (|6.ip shows that the action of G on M induces another group action of G on 
E n . This yields a partition of E n into orbits which can be described by the original 
action of G on M. For a given element (Ai(X)) in E n , we set 

Tij — AijAj/Ai. 



Theorem 8. Suppose that the coefficients Ai, for n — 1, . . . ,n in equation (|1.2p 
are non vanishing and that exactly p of the functions A^j in the expression of the 
generators V^i 1 of the first prolongation of G in (|4. 3[) vanish identically. Then two 
differential equations (Aj(X)) and (Bi(X)) are equivalent, i.e. they belong to the 
same orbit of the first prolongation of G if and only if 

T$(X) = 7${X) 
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for all of the n(n — 1) — p nonzero such functions TA and . 

Proof. If the coefficients Ai of (|1.2[) are non vanishing, then diag {Ai, . . . , A n } has 
constant rank n, and by the expression of the generators in (|4.3[) , G acts semi- 
regularly. More over, the expression of the corresponding invariant functions in 
(|4.5[) shows that this action is regular. Since the invariants of G are actually 
the invariants of the induced group action of G on E n , the result follows from 
a theorem of [TOl Theorem 2.34] stating that for regular group actions, two points 
lie in the same orbit if and only if they take on the same values under all invariant 
functions. □ 

7. Concluding Remarks 

It follows from Theorem [8] above that all differential equations (Ai(X)) where 
Ai = Ai(x % ) depends on x % alone are equivalent. The expression of the generators 
in (|4.3p shows that the action of G restricted to this family of equations has no 
invariant of any order. 
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